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Abstract 



In this paper, we study the Arneodo-Coullet-Tresser map F{x, y, z) = {ax — 
b{y — z),bx + a{y — z),cx — dx k + ez) where a, b, c, d, e are real with bd and 
k > 1 is an integer. We obtain stability regions for fixed points of F and symmetric 
period-2 points while c and e vary as parameters. Varying a and e as parameters, 
we show that there is a hyperbolic invariant set on which F is conjugate to the full 
shift on two or three symbols. We also show that chaotic behaviors of F while c 
' and d vary as parameters and F is near an anti-integrable limit. Some numerical 

results indicates F has Hopf bifurcation, strange attractors, and nested structure of 
invariant tori. 

1 Introduction 

In this paper, we consider discrete dynamical systems induced by maps F : M 3 — ► M 3 of 
the form 

F{x, y, z) = {ax — b{y — z),bx + a{y — z),cx — dx k + ez), (1) 

where a, b, c, d, e are real parameters with bd ^ and k > 1 is an integer. These maps 
are introduced by Arneodo, Coullet and Tresser being motivated by the study of strange 
attractors in a family of differential equations on M 3 with homoclinic points of Shilnikov 
type, refer to [5]. Their numerical computations showed some interesting phenomena in 
dynamical behavior of these maps (unpublished); in particular, for the one-parameter 
family of maps (pQ) with a = 0.6, b = 0.5, d = e = l,k = 3 and c as the parameter, they 
discovered (what may be regarded as) period doubling cascade and a strange attractor. 
Rigorous results on characterization of periodic points of (OQ), as well as existence of 
topological horseshoes for a region of parameters were obtained in plj. Nevertheless, since 
dynamical behavior of maps (CD) is rich and multifarious in several regions of parameters, 
there remains many interesting problems in rigorous understanding the dynamics and 
describing the bifurcation structure of ([T]); some of these problems will be considered in 
the present paper. 



"Institute of Mathematics, Academia Sinica, Taipei 115, TAIWAN, dubs@math.sinica.edu.tw 
^Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, TAIWAN, 
mcli@math.nctu.edu.tw 

■'■Department of Mathematics, Nizhny Novgorod State University, Nizhny Novgorod, RUSSIA, 
malkin@uic . nnov . ru 



1 



We will call the maps ([I]) the Arneodo-Coullet-Tresser maps, or ACT maps for short. 
Note that the Jacobian determinant of an ACT map F is a constant; namely, I d Q{xyz) I = 
(a 2 + b 2 )e. Furthermore, if e 7^ then the map F is a diffeomorphism with the inverse 

—bx + ay 



F (x, y, z) = [x, + z, z) 

a z + tr 



where x = — — — and z = . So, t is a polynomial automorphism of 

cr + b 1 e 

M 3 and thus, in the case of ACT maps, the Jacobian conjecture holds true (from this 
point of view, the ACT maps can be compared with generalized Henon maps, which are 
also polynomial automorphisms, see [6j [11] about the history and results concerning the 
Jacobian conjecture, see [3| H2])- 

When studying dynamical systems on noncompact manifolds, it is desirable to be 
sure that the nonwandering set of the system is compact, in which case one might restrict 
oneself to this nonwandering set in order to consider nontrivial dynamical behavior. In [7], 
a sufficient condition for the nonwandering set of polynomial maps to be compact was 
obtained and an estimate for the size of the box containing the nonwandering set was 
given. These results being applied to the case of ACT maps imply (see Proposition 
15]) that the nonwandering set £l(F) (as well as the set of bounded orbits) of the ACT 
map F lies in the box 

x, y, z el : x < M, y < — ^ M, z < — —M 

, j \b\ " 1 - \b\ 

where 

M= k ~ 



1 \a 2 e + b 2 e\ + \a 2 + b 2 - bc + 2ae\ + \2a + e| + 1 



\bd\ 

Note that the above result of the nonwandering set is a generalization of a result in [5] 
which asserts that the set of periodic points of the ACT map is contained in the box. 

We will discuss (see Sections 3 and 4) chaotic dynamical behavior of the ACT maps 
in some parameter region. Usually chaotic behavior of a dynamical system is associated 
with positive topological entropy. For maps not uniformly continuous with non-compact 
domain one needs to be very specific to define the topological entropy h top {F) to be 
htop{F\n(F)) ■ On the other hand, since F is a polynomial automorphism of K 3 , it can 
be easily extended to the homeomorphism F of the canonical one point compactification 
R 3 := R 3 lJ{oo} with F(oo) = 00 and so one can also define ht op {F) to be ht op {F) and 
the two definitions of the topological entropy agree because h top {F) = h top {F\^p^) and 
00 is an isolated nonwandering point for F. 

The paper is organized as follows. In Section 2, we obtain stability regions for fixed 
points of the ACT maps and symmetric period-2 points for two-parameter family of the 
ACT maps with c and e as the parameters. These stability regions are used later in Section 
5 to explain period doubling bifurcation and discrete Andronov-Hopf bifurcations. We also 
observe that the intersection of these stability regions is a point at which the Jacobiant 
determinant of the ACT map takes value 1 (and so the system becomes conservative) and 
which gives rise to a nested family of invariant two-dimensional tori, and as numerical 
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study shows, this nested structure of invariant tori persist for some interval of parameter c 
with parameter e fixed (so that that system remains conservative). In Section 3, we study 
chaotic behavior of the ACT maps for two-parameter family with a and e as parameters. 
We prove that for a and e small there is a hyperbolic invariant set which is conjugate to 
the full shift on two or three symbols depending on the evenness of the number k in ([1]). 
In Section 4, based on a result in [8], we study the chaotic dynamics of the ACT maps for 
two-parameter family with c and d as parameters. While c — > oo and d —>■ oo in such a 
way that - = constant > 0. Let A = -. Then for all |A| sufficiently small, the ACT map 
F — Fx has a closed invariant set A^ such that F\\Ax is conjugate to the full shift on either 
two or three symbols depending on whether k is even or odd, respectively. In Section 5, 
we discuss numerical results on bifurcations of the ACT maps including Hopf bifurcations 
and resonances, appearance and changing of invariant circles, strange attractors, nested 
structure of invariant tori. In the Appendix, we derive the stability criterion for 3x3 
real matrices, which is used in Section 2 (and is more convenient for our purposes than 
the Shur-Cohn criterion). 

While studying parameter-dependence property of F, we will denote F as F c (resp. 
F e ) to stress that c varies (resp. e varies) as other parameters are fixed. Similarly, we 
denote F c e when c and e are the only parameters that vary. 



2 Stability regions for fixed points and symmetric 
period-two points 

It is clear that the origin is a fixed point of F. By solving F(x, y, z) = (x, y, z), we obtain 
that for even k, the map F has a unique nontrivial fixed point at 

a 2 + b 2 -a (a - l) 2 + b 2 , 
Pi = {xi, £ xi, - Xx) 

where X\ = k ~\J / fec ~( 1 ~ e )K^~ 1 ) anc j f or case when k is odd and fec ~( 1 ~ e )K°~ 1 ) + b 1 > 
0, the map F has exactly two nontrivial fixed points at dbpi. Moreover, by solving the 
system of equations F(x, y, z) = (—x, —y, —z) and F(—x, —y, —z) = (x, y, z), we get that 
F has periodic points of period 2 symmetric to the origin, say ±p 2 , if and only if k is odd 
and 6c ~( 1+e )[fc+ 1 ) +b 1 > 0; under these conditions p 2 is given by 

-a 2 -6 2 -a -(a + lf-b 2 . 
P2 = {x 2 , 7 x 2 , t x 2 ), 



with x 2 = k y C m • We call ±p 2 (when exist) the symmetric period-two points 

of F. 

In next section, we will be concerned with stability regions for fixed points and period- 
2 points of the ACT maps, i.e., regions in the parameter space for which these points are 
stable. More precisely, we use the following definitions. 

Definition 1. Let x i— > G v (x),x G X C M m , v G V C M! be a family of C 1 maps with 
parameter v, and let x* be a fixed point of G v for each v G V . A subset J C V is called 
the stability region for the family of fixed points x* if for each v G J, G v is linearly 
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stable at x* and for each v ^ J, G v is not linearly stable at x* ; i.e., for each v G J all 
eigenvalues of the Jacobian matrix ^G v (x*) lie in {z G C : \z\ < 1} and for each v ^ J 
there is an eigenvalue of the Jacobian matrix ^G v (x*) in {z G C : |z| > 1}. 

Similarly, we can define a stability region for a family of a periodic orbit of period n 
by replacing G v by G™. 

For the ACT family F, the Jacobian matrix of F at a point (x, y, z) is 



dF(x,y,z) 
d(x,y,z) 



a —bb 
b a —a 



c — kdx k 1 e 
and its characteristic polynomial is 

P(A) = A 3 - (2a + e)A 2 + [a 2 + b 2 + 2ae - be + kbdx^X - (a 2 + b 2 )e. 

d F (x , t/, z\ 

Note that the determinant of the Jacobian matrix of F is constant, that is, I — — ; — - — ^— I - 

d{x,y,z) 

{a 2 + 6 2 )e, and if e 7^ then the map F : M 3 — > M 3 is a diffeomorphism with the inverse 

p-i/ \ /- -bx + ay A 
^ (g,y,^) = (a, 2 , , 2 

+ 

ax + by z — cx + dx k „ 
where x = — — — and 2; = . So, F is a polynomial automorphism or K 

a 2 + o 2 e 

and thus, in the case of ACT maps, the Jacobian conjecture holds true (about the history 
and results concerning the Jacobian conjecture, see [3], [T2]). 

To determine the local stability of the fixed points, we need to know whether eigenval- 
ues of the Jacobian matrix of ACT map lies inside the unit circle in the complex plane, 
i.e., whether roots of the characteristic polynomial of this matrix lies inside the unit circle. 
We will use the following definition. 

Definition 2. We say that a polynomial P over M is stable if all of its roots lie in 
{z G C : \z\ < 1}. 

Note that polynomials of this type of stability are called sometimes discrete stable or 
Shur stable to distinguish from the situation, when all roots of polynomial have negative 
real parts (this applies usually to characteristic polynomials of linearization of vector 
fields at stable fixed points), in which case the term Hurwitz stability is sometimes used. 
There are several criteria to determine stability of a polynomial in terms of its coefficients 
(e.g., the Shur-Cohn criterion, see [TO]). For our purposes (taking in mind bifurcation 
parameters of ACT maps) it will be convenient to use the following one, which is derived 
in the Appendix. 

Proposition 3. A polynomial P(X) = A 3 + AX 2 + BX + D G M[X) is stable if and only if 

\D\ < 1 and max{-P(l), P(-l)} < < a, 
where a = &(P) := —D 2 + AD-B + 1. 
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We apply now the above criterion to a one-parameter family of polynomials, which 
are the characteristic polynomials of the Jacobian matrix of ACT maps. The proof of the 
following proposition is postponed to the Appendix. 

Proposition 4. Let Pp(X) = A 3 + AX 2 + (B + (3)X + D be a one-parameter family of 
polynomials in M[A] with parameter /3 G R. Then 

1. The following three statements are equivalent: 

(a) \D\ < 1 and \ A — D\ < 2. 

(b) There exists a (3 such that all roots of Pp lie in {z G C : \z\ < 1}. 

(c) There is a unique open interval I cl, which will be called the stable interval 
of the family Pp, such that for each (3 e I, Pp is stable and for each (3^1 
there is a root of Pp in {z G C : \z\ > 1}. 

2. The stable interval of Pp {when exists, see item 1) is equal to (a,a ), where a = 
max{— P (l), P (— 1)} and a = —D 2 + AD — B + 1. Moreover, for (3 = a the 
polynomial Pp has two complex conjugate roots in {z G C : \z\ = 1}\{— 1, 1}. 

Now we are in position to study the stability regions for the map F. In the following 
Theorems [M3 we deal with stability regions J tr (-fc) and Jtr(-^c,e) (recall that the subscripts 
here indicate the only parameters that vary). Note that if we find some functions /i(e), 
/2(e) of variables e, such that Jtr(Fc,e) = {( e ) c ) £ ^ 2 : /i( e ) < c < /2(e)} then we will 
have for any c with J tr (P c ) 7^ 0, that J tT (P c ) = {c£l:/i(e) < c < f 2 (e)} , i.e., Jtr(-^c) is 
described by the same inequalities. So in this case it is enough to give the corresponding 
formulas for Jtr(F c ,e) only. 

Theorem 5 (stability regions for the trivial fixed point). Let F be the ACT family with 
b 7^ 0. Let J tr (F c ) (resp. Jt r {Fc,e)) denote the stability region of the origin for F c (resp. 
forF c ^ e ). Then 

1. J tr (F c ) ^ if and only if 

-1 < (a 2 + b 2 )e < 1 and 2a - 2 < (a 2 + b 2 - l)e < 2a + 2. (2) 

2. For F c ^ e , the following two statements hold: 

(a) Suppose a 2 + b 2 — 1 < 0, then Jt r (F Cje ) 7^ 0. 

(b) Suppose a 2 + b 2 — 1 > , then Jt r (F c ,e) 7^ if and only if 

f 2a -2 -2a - 2 
max { a^ + b 2 -! 1 a 2 + b 2 -l 

3. If J tr {F C) e) ^ 0, then 

Jtr(F c ,e) = {(e, c) G M 2 : max{-ci(e), c_i(e)} < -6c < c(e)}, 

where 

Cl (e) = (1 - e)[{a - l) 2 + b 2 ], c. x {e) = -(1 + e)[{a + l) 2 + b 2 ] 
and c(e) = -(a 2 + b 2 - l)[(ae - l) 2 + 6 2 e 2 ]. 

(See Figure 1 and the remarks below it. ) 



< 



+ b 2 ' 



(3) 
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Proof. Since b ^ 0, we may use (3 = —be as parameter. Then the characteristic polynomial 
of the Jacobian matrix of F at the origin is 

p p (\) = A 3 - (2a + e)A 2 + (a 2 + b 2 + 2ae + (3)X - (a 2 + b 2 )e. 

By applying item 1 of Proposition H] to the above family Pp with parameter [3, we obtain 

that equation (4) is equivalent to that the stable interval of Pp is not empty and so is 

equivalent to that the stability region of F is not empty. Item 1 is completed. 

The inequalities a 2 + b 2 — 1 < and b ^ imply \a\ < 1, so the numbers 2a — 2 and 

2a + 2 are of opposite signs, and therefore one gets item 2(a) from item 1 by taking e = 0. 

• „/, n , , ±1 2a -2 , 2a + 2 

For item 2(b), let us denote e± = - ■ — , e x = - and e r = - -. 

a 1 + tr a 1 + (r — 1 a z + b z — \ 

The existence of e satisfying condition (4) is equivalent to the fact that the two intervals 
(e_,e + ) and (e/,e r ) overlap, i.e., e; < e + and e_ < e r . It is easy to see that the last two 
inequalities are the same as (5). 

The stable interval of Pp is given by max{— Pq{1), Pq(— 1)} < j3 < do, (for the defini- 
tion of do see item 2 of Proposition HJ). By evaluating the values Pq(1), Pq{— 1) and do 
and using the fact that (3 = —be, item 3 is proved. □ 

Next we study the stability region of the nontrivial fixed point (s), say pi for even k and 
±Pi for odd k. Let X\ be the x-coordinate of pi, then one easily gets that %\ = ' 



bd ' 

where c x = (l-e)[(a-l) 2 +6 2 ] (the same as in item 3 of Theorem [5]). Similar to Theorem [5], 
applying Proposition S] to the characteristic polynomial of the Jacobian matrix of F at 
the nontrivial fixed point (s) p\ (or ±pi), which is 

Pp{\) = A 3 - (2a + e)A 2 + (a 2 + b 2 + 2ae + /3)\ - (a 2 + 6 2 )e, 

where (3 = —be + kbdx^ 1 , we have the following result. 

Theorem 6 (stability regions for the nontrivial fixed point (s)). Let J n on{F c ) {resp. 
J non {F c>e )) denote the stability region of the nontrivial fixed point(s) for the map F c (resp. 
for {F c>e )). Then 

1. If k is odd and bd < then J n0 n(F c ) = J n0 n(F c ,e) = 0- If k is even or bd > 0, 
then Jnon(F c ) ^ if and only if Jt r (F c ) ^ 0; moreover J n0 n(F Cte ) ^ if and only if 

Jtr{Fc,e) ^ 0- 

2. If Jnon{F c , e ) ^ 0, then 

Jnon{F c , e ) ={(e,c) G M 2 : 

fcci(e) + c(e) . j. , . fec 1 (e)+c_i(e) 1 i 
^— j < -6c < mm{- Cl (e), ^— }|, 

where c(e), Ci(e) and c_i(e) are m item 3 of Theorem^ 

3. For F Cte! the point (— ci(e),e) zs on i/ie boundary of the stability region of the non- 
trivial fixed point(s) if and only if it is on the boundary of the stability region of the 
origin. 

See Figure 2 for the regions of J non (F ce ) aside with J tr (F ce ). 
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For the symmetric period-two points at ±p 2 of the ACT family F with odd k, we have 

DF 2 (p 2 ) = DF(p 2 ) ■ DF(-p 2 ) = \DF{p 2 )} 2 

and so the eigenvalues of DF 2 {p 2 ) are the square of the eigenvalues of DF{p 2 ). It fol- 
lows that the stability of ±p 2 is determined by whether the characteristic polynomial 
of DF(p 2 ) has all roots in the open unit disk. Let x 2 be the x-coordinate of p 2 , then 

x<1 = k ~\J bC+ bd~^ where c_i is in item 3 of Theorem O Applying Proposition @] to the 
characteristic polynomial of DF(p 2 ), which is 

P p {\) = A 3 - (2a + e)\ 2 + (a 2 + b 2 + 2ae + (3)X - (a 2 + 6 2 )e, 

where (3 = —be + kbdx 2 _1 , we have the following result. 

Theorem 7 (stability regions for the symmetric period-two points). Let the ACT family 
F have the symmetric period-two points ±p 2 , i.e., k is odd and bdibc + c_i(e)) > ; and 
let J sym (F c ) (resp. J sy m{F Cfi )) denote the stability region of the symmetric period-two orbit 
for F c (resp. for F C)e ). Then 

1. Ifbd<0 then J sym (F c ) = J sym (F C:e ) = 0. Ifbd > 0, then J sym (F c ) ^ if and only if 
Jtr{F c ) ^ 0/ moreover J non {F^ e ) ^ if and only if J ir (F Cje ) ^ 0. 

2. If J sym {F C) e) ^ 0, then 



fcc_i(e) - c(e) . /cc_i(e) + ci(e) 
— < -be < mm{ — , c_i(e)} 

where c(e), Ci(e) and c_i(e) &e m ztem 3 of Theorem^ 

3. For F c ^ e , the point (c_i(e),e) is on the boundary of the stability region of the sym- 
metric period-two points if and only if it is on the boundary of the stability region 
of the origin. 

See Figure 3 for the regions of J sym (i 7 'c,e) together with J tT (F c>e ) and J n on(F c>e ). 

Remark 8. Note that for k odd, the condition bc ~( 1+e )(^ +1 ) +b ) > q j or existence of 
symmetric period-two points is satisfied precisely in the half-plane of (e, c) -plane which 
is boundary by the straight line c = c_i(e) and lies to another side from this line than 
Jtr(F c , e ) the stability region of the trivial fixed point. So by using item 3 of Theorem 7 one 
map associate appearance of these period-two points with period doubling bifurcation of 
the trivial fixed point, which occur at the common segment do(Jt r (F Cje )) fl clo(J sym (F C)e )) 
of the line c = c_i(e) (recall that one of the eigenvalues of Jacobian matrix at the trivial 
fixed point is —1 as far as parameter (e, c) belong to this line). 
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3 Hyperbolic structure and infinitely many periodic 
points 

In [5], the first author of this paper proved that the set of periodic points of the ACT 
family is contained in a bounded box. In [7], the other two authors extent the result as 
follows. 



Proposition 9. The nonwandering set of the ACT family F lies in the bounded box 

+ 2|a| + l , a 2 + b 2 



s , , , , « 2 + ^ 2 + 2|a| + l_ , a 2 + b 2 +\a\ nr 
x, y, z) e K : \x\ < M, \y\ < -— ^ M, \z\ < — —M 



where 

_J\a 2 e + b 2 e\ + \a 2 + b 2 - bc + 2ae\ + \2a + e\ + 1 



M 



\bd\ 



Remark 10. Note that if f is a homeomorphism from M > onto itself then f can be 
easily extended to a homeomorphism f of the canonical one point compactification M. m : = 
M > |J{oo} with /(oo) = oo. Therefore the topological entropy h top (f) can be defined (see 
ffijll) as the usual topological entropy of a homeomorphism of a compact set. On the other 
hand, using the results of JW on boundedness (and hence compactness) on fl(f) under our 
assumptions, one may define h top (f) to be h top (f\ n{f) ) . Since h top (f) = h to p(f\ n ^) and 
oo is an isolated nonwandering point for f , the two definitions of the topological entropy 
of f agree. 

The following theorem shows the existence of horseshoe. 

Theorem 11. Let F afi be the ACT family with parameters a and e. Suppose b,c,d,k 
satisfy 

(i) \b\ > 1 and \b - c\ > 1 
and one of the following: 
kb 2 



(iil) k is even, be > 



k-V 



(U2) k is even, 

rC 1 

kb 2 _ 

(m3) k is odd, bd > 0, be > (1 + vk), 

k — 1 

kb 2 _!/— 
(UA) k is odd, bd < 0, be < -(1 - k vk). 

Then for each pair of parameters a and e which are sufficiently close to 0, F a<e has an 
invariant set, say A a>e , such that the following properties hold: 

1. F ae has a hyperbolic structure on A ae . 
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2. A aie is a Cantor set. 



3. -F aie |A aie is topologically conjugate to the two-sided shift on two (resp. on three) 
symbols for even k (resp. for odd k). 

Proof. By the persistence of hyperbolic invariant sets, we only need to consider the cease 
when a = e = 0. We divide the ACT family into 8 cases depending on the signs of d and 
be for even k and on the signs of b and d for odd k. Since the proofs of all these 8 cases 
are similar, we only give a proof of the case when k > 1 is even, d > 0, and be > 0. Since 

F 0i0 : (x, y, z) ^ (x, y, z) = (-b(y - z), bx, cx - dx k ), 

by eliminating x among the expressions of y and z, we obtains that the image of R 3 under 
F 0i0 is the two-dimensional surface 

r c d 

T = Ux,y,z) : x ER,z = -y - —y k 

dz c d 

By letting — = for z = -y — —ry k , we have that the surface T has the line 

J 6 dy b b k 



(x,y*,z*):xeR,y*= ^—^ = ^- ^ — 

as the set of extreme points. By letting z = 0, we get that T intersects the (x, y)-plane at 
the lines y = and y — r* where = y — - — . It follows from the assumptions d > 

and 6c > that r* > > and > 0. Also the assumption be > implies z* > y*. 

K — 1 

Since x = — b(y — z), F o maps the infinite strip 

S = {{x, y, z) : x e M, -r* < 2 - y < 0} 

onto the surface 

r c d 

Foo(S) = < (x, y , z) : z = -y — j-ry k , x lies between and — br* 
y b b k 

which is a subsurface of T. Therefore, S D F 0i o(5 l ) consists of two disjoint connected sets 
whose projection to the (y, z)-plane are shown in (i) of Figure 4; (ii) of Figure 4 is for the 
case when k is odd. 

Next, we study Fqq(S) H S, where Fq'q 1 denotes the preimage of S under F 0; o- Con- 
sidering z — y = — and z — y = 0, one has the equations dx k + (b — c)x — r* = 
and dx k + {b — c)x = 0. The two equations give four distinct real roots in x, namely 
Ai < A 2 < A 3 < A 4 . Then 

^oToOS) = {(x,y,z) : A x < x < \ 2 ,y e R,z G R} 

\J{(x, y,z):\ 3 <x<\ i ,yeR,ze R}. 



9 



Let A = Fq Q (S). We need to check the hyperbolicity conditions for points in 

a small neighborhood U of Fqq(S) PI S. First of all, we make a change of coordinate via 
u = y — z and v = y + z. Then in the new coordinate (x, u, v), we have 



F 0)0 (x,u,v) 



-bu, bx — cx + <ix fc , 6x + cx — dx k ). 



Let x = (x,y) G R 2 and write F 0i0 (x, m, = (G(x, v), H(x, v)), where G : R 3 — > R 2 and 
H : R 3 — > R. Then the sufficient conditions of hyperbolicity are the following (refer to 

CD: 



(iH <;L > 



dv 



< 1, 



i-IKS-II 

(i-ii(i)-i 



9/7 
(1- 



> 2 



ii— r— v l \\ 

9x ^ 9x ' 
,,dH ,dG.-i t 

> hr-(ir) 

ax v ax ' 



dG 

dv 



(5) 
(6) 

(7) 



where || ■ || = sup | ■ |. 

(x,v)GC7 

The Jacobian matrix of F in the new coordinate is 



'dG 


dG' 




9x 


dv 




dH 


dH 




. <9x 


dv . 





-b 

b-c + kdx k - x 

b + c- kdx k ~ l 



Since 

dv 



|^|| = 0, (7) and (8) are clear and (9) follows from (6). We only need to 



check (6). By the definition of Aj's and the derivative of f(x) = dx k + {b — c)x, we have 
that 

min{|6 - c + kdx k ~ l \ : x G [A 1; A 2 ] U [A 3 , A 4 ]} = \b — c\. 
Thus (6) follows from the assumptions \b\ > 1 and \b — c\ > 1. □ 



4 Chaos near an ant i-integr able limit 

In [8], the second and third authors consider solutions of families of difference equations 
^x{Vni Vn+i, ■ ■ ■ , Vn+m) = 0, u G Z, with m+1 real variables for parameters A near those 
exceptional values Ao for which the difference equation function depends on only one 
variable: $>\ (yo, . . . ,y m ) = <p(yN) with some integer < N < m and function if. It is 
proved there that if (p has k simple zeros, then among solutions for parameter values close 
to exceptional ones , there are "topological k- horseshoes" , i.e., solutions for which the 
restriction of the shift map is conjugate to the two-sided full shift on k symbols, and if the 
difference equations correspond to smooth diffeomorphisms on R > , then the topological 
£;-horseshoes persist among orbits of such diffeomorphisms. For the exceptional parameter 
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values, no real maps are correspondent, and such a situation reminds the so called anti- 
integrable limit approach inspired by Aubry and Abramovici's anti-integrable concept 
in 0; see also 01 El El. 

For an initial point p = (x Q ,y ,z ), denote the n-th iteration of p under the ACT 
map F by (x n , y n , z n ). Then for any nGZ, we have the following system consisting of 7 
equations (the index i for labeling the equations bellow takes the values 1 and 2): 

X n +i = ax n+i-l ~ ^{Vn+i-l ~ z n+i-l)j (1) V 

Un+i = ^> x n+i-l + a (lln+i-l ~ z n+i—l)i (2, i) 

z n+i CX n -\-i—i dx n ^i__-[ ~\- C2Vi+j— 1; (3,i) 

x n+3 = ax n+2 - b(y n+2 - z n+2 )- (1, 3) 



From equations (1,1) and (2,1), we can express y n+ \ in terms of x n and x n+ \. Then 
plugging y n+1 into equation (1,2), we can express z n+ i in terms of x n , x n+ i and x n+2 . 
The expressions are 



a 2 + b 2 a 
~b~ : 

a 2 + b 2 2a 



Vn+l , X n i -Kn+l-i 

b b 





Similarly, from equations (1,2), (2,2) and (1,3), we can express z n+2 in terms of x n+ i, 

x n+2 and x n+3 . Then, plugging the expressions of z n+ \ and z n+2 into equation (3,2), we 
get the difference equation 

k a 2 e + b 2 e a 2 + 6 2 — be + 2ae 2a + e 1 

dx n +i H j Xn H ; ^n+l ; + = 0. 





Suppose the parameters a, 6 and e to be fixed and let c — > oo, d — > oo in such a way 
that - = constant := A > 0. Denote A = -, then d = 4 and for the ^-coordinate of orbit 

c c ' A 

under the ACT map / = j\, we have the difference equation 

fc _ x / a 2 e + b 2 e a 2 + b 2 + 2ae 2a + e 1 \ _ n 

^n+u^n+i — A) + A I a; n H x n+ i — x n+2 + -^%n+3 J — 0. 

(9) 

The difference equation ([9]) corresponds to the map fx for A ^ because between solutions 
x = (a; n )^ = _ 00 of ([9]) and full orbits p = (p™)^.^ of /a, we have a conjugacy inp, 
given by 

a 2 + 6 2 a a 2 + 6 2 2a 1 

Pn \Xni , X n —\ —X ni - X n —\ , X n -\- ~X n j^2j i 

b o o oo 

while the inverse p i— > x is given by x n = 7Ti(p n ). For the limit value Ao = of parameter 
(in this case, it is an anti-integrable limit; refer to [2]), we have 

if)(x n+ i) := x n +i{x*+\ - A). 

Thus the function ip has at least two simple zeros; more precisely ip has two simple zeros 
{0, k ~\fA} when k is even, and three simple zeros {0, ± k ~\/~A} when k is odd. By applying 
Theorem 3 of [S], we have the following results. 
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Proposition 12. Suppose that for the family of the ACT maps ([IT) , the parameters a, b, 
and e are fixed while c — > oo and d — > oo in such a way that - = constant > 0. Let A = -. 
Then for all |A| sufficiently small, the ACT map F = F\ has a closed invariant set A\ 
such that F\\A\ is conjugate to the full shift on either two or three symbols depending on 
whether k is even or odd, respectively. 

There are also other parameter routes in families of Arneodo-Coullet-Tresser maps, 
for which the above arguments apply. Namely, if we consider a and c to be fixed while 
e — > 0, and b — > oo and d — > oo in such a way that 4 = constant > 0, then similarly to 
the lines of the previous proposition, we will have the following. 

Proposition 13. Suppose that for a family of the ACT maps ([!]), the parameters a and 
c are fixed while e — > 0, and b — > oo and d — > oo in such a way that 4 = constant > 0. 
Let X = ~. Then for all sufficiently small \e\ and \ X\, the ACT map F = F e ^\ has a closed 
invariant set A EjA such that F ejA |A eiA is conjugate to the full shift on either two or three 
symbols depending on whether k is even or odd, respectively. 



5 Numerical results 

Fix a, b, d, k such that the stability regions of the origin, the nontrivial fixed point p\ 
and the symmetric period-2 point p2 are not empty. For any e we denote the boundaries 
of the stability regions as follows: c tr (e) = c_i )tr (e) = £ ^zj^, c sym (e) = hc ~_^Zx)^ , 

C-l,non(e) = ^^fct-'i") 1 ^ ; and C non (e) = ^if^-l)^ ' wllere K e )i C l( e ) aIld C -l( e ) are in item 

3 of Theorem 5; see Figure 5. Let e± = =fc a2 ^_ fc2 , then c tr (e ± ) = — — — 1 ^^ a ^ ± ^ +b e± ^ 
and the Jacobian matrix of F Cl t e± \ e± at all fixed points has all eigenvalues lying on the 
unit circle. For e with e_ < e < e + , the Jacobian matrix of F^r e \ e (resp. F^ ( e ), e 
and -PgnonCe)^) at the trivial fixed point (resp. at the symmetric period-2 point and at 
the nontrivial fixed points) has simple complex eigenvalues lying on the unit circle and 
the Jacobian matrix of F c _ ltI < e \ e (resp. F c _ 1 non ( e ), e ) at the trivial (resp. nontrivial) fixed 
points has simple real eigenvalue at —1. Numerical results indicate the following: 

1. The family F Cje undergoes Hopf bifurcations at (c, e) equal to (c tr (e), e), (c sym (e), e) 
and (c non (e),e) for almost all e with e_ < e < e+. As the value of c decreases, 
the former bifurcation is subcritical and the invariant circle (before disappeared) is 
saddle type while the latter two are supercritical and the appeared invariant circle is 
asymptotically stable. It can be shown that resonance cases take places for several 
parameters; e.g., for c equal to c tr (e) and c non (e), F Cj6 has resonance (1 : 4) at 
e = a 2+b2_i and resonance (1 : 3) at e = a 2^ 1 _ 1 - See Figure 5. 

2. The family F c>e undergoes period-doubling bifurcations at (c, e) equal to (c_i jtr (e), e) 
and (c_i jnon (e), e) for all e with e_ < e < e, where e with c_i >tr (e) = c_i >non (e). In 
fact, one can show that the symmetric period-two points ±p2 bifurcates from the 
origin at (c_i )tr (e), e). See Figure 5. 

3. If e = e±, then for each c near c tr (e±) with d(c — c tr (e±)) > 0, the map F c has 
infinitely many invariant tori (or double tori); one lies inside another and none of 
these tori is asymptotically stable. See Figure 6(i) with a = 0.2, b = — 1.4, c = 
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— 0.94, d = — l,e = 0.5 and k = 3 and Figure 6(ii) with a = 0.6, b = 0.8, c = 
— 0.01, d = — l,e = —1 and k — 4. For the former case, one invariant circle is 
drawn and it indicates that all other invariant tori are bifurcated from the invariant 
circle through a degenerate Hopf bifurcation by taking a Poincare section. For the 
latter case, each double torus consists of two tori which are mapped to each other 
alternatively and is invariant under the map F c . 

4. For some parameters, e.g., a = b = 0.5, d = e = 1, and k = 3, the family F c under- 
goes sequences of wind-doubling bifurcations of invariant circle, and has coexistence 
of stable invariant circles and strange attractors for some parameters c. See Figure 
7. 

5. For some parameters, e.g., a = e = 0.01, b = 1.1, c = 3.6578, d = 1, and k = 3, the 
family F has a strange attractor. See Figure 8. 



Appendix: Proofs of Propositions [3] and [4 

Proof of Proposition^ First, we consider the case when D = 0. Let Q(X) = A 2 +AX+B, 
then -P(A) = XQ(X). For the "only if part of the proposition, the stability of P implies 
that the product of two roots of Q has absolute value less than one and so B < 1, i.e., 
a > 0. Together with the fact that the leading coefficient of Q is positive, we have that 
P(l) = Q(i) > and -P(-l) = Q(-l) > 0. For the "if part, by the assumptions, we 
have that B < 1, 1 + A + B = Q(l) = P(l) > 0, and l — A + B = Q(-l) = -P(-l) > 0. 
Adding the last two inequalities gives that B > — 1 and so together with the first inequality 
we obtain \B\ < 1, that is, the product of two roots of Q has absolute value less than one. 
It follows that the two roots of Q lie inside the unit circle because of the following facts: 
Q(l) = -P(l) > 0, Q(— 1) = — P{— 1) > 0, and that the leading coefficient of Q is positive. 
The proof of the proposition in the case when D = is completed. 

Next, we consider the case when D^0. Before the proof, we give some preparations. 
Define P a (X) = P{X) + aX as a one-parameter family of polynomials with parameter 
a G R. Then for any fixed Ao ^ 0, P a (Ao) has the monotonicity property in a: if 
A > (Ao < 0, A = respectively), then P a (Ao) increases (decreases, equals to D, 
respectively) as a increases. It is clear that there is a unique a, namely a = — P(l), 
such that P a (l) = 0. Similarly, P a (— 1) = for a unique a, namely a = P(— 1). Since 
the product of the three roots of P a is —D, it follows that there is a unique a, namely 
a = P ^~ D - > = a, such that P a has two roots whose product is equal to one. In particular, 
if P a has two complex conjugate roots on the unit circle, then a = a. For convenience, 
we denote a\ = — P(l), a_i = P(— 1) and a max = max{ai, a_i, &}. Note that the 
numbers a\, ct_i, a partition the real line into disjoint open intervals and for any a, a' in 
the same interval of this partition, we have that P a and P a > are either stable or unstable 
simultaneously. Indeed, this follows from the fact that roots of polynomials depends 
continuously on the coefficients. Because of this, we call ai, a_i and a the bifurcation 
values of the family P a . 

Now, we prove the "only if part. The first statement \D\ < 1 follows immediately 
from the stability of P. To prove the second statement, we need the following two claims: 
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Claim (i) P a is unstable for all a > a max . Suppose the claim is not true. Then there 
is a' > a max such that P a > is stable. Thus for each a > a max the polynomial P a is stable 
and so X a := min{A G (—1, 1) : A is a root of P a } is well defined. Thus 

|aA Q | = | - P(X a )\ = \X 3 a + AX 2 a + BX a + D\ < 1 + \A\ + \B\ + \D\. 

It follows that X a — > as a — > +oo. Thus for sufficiently large a, the product of the 
two roots of P a other than X a is whose absolute value is large. This contradicts the 
stability of P a . 

Claim (ii) P a is unstable for all a < max{ai,a_i}. Indeed, by taking Ao = 1, the 
monotonicity property in a of P a (A ) implies that P Q (1) < P ai (l) = for all a < a.\. On 
the other hand, the positivity of the leading coefficient of P a implies that P a (X) > for 
some A large. Thus for all a < ai the polynomial P a has a real root bigger than 1 and 
so P a is unstable. The proof of the fact that P a is unstable for all a < ct_i is similar by 
taking Ao = — 1. 

Combining the claims (i) and (ii), we get that if a < maxjai, then P a is unstable 
for all aGi In other words, if P a is stable then max{«i, < a < a. In particular, 
by the stability of P a with a = 0, we have that max{«i, <x_i} < < a. 

Next, we prove the "if part. Since the only bifurcation values of P a with parameter a 
are ai, a_i and a, it is sufficient to show that P a is stable for some a' with max{«i, a^i} < 
a' < a. By the assumption \D\ < 1 and the definition of a, for a = at the polynomial P a 
has a real root at — D which lies in (—1, 1). Since a > max{«i, a_i} and the monotonicity 
property of P a (Ao) with Ao = 1 and —1 imply P&(— 1) < < Pq(1). Now take an a' 
slightly less than a so that max{«i,a_i} < a' < a and P a r(—1) < < P a >(l). By the 
monotonicity property of P a {Xo) with Ao = —D, we get that P a > has a real root, say A', 
with | — D\ < |A'| < 1. Let Ai, A2 denote the other two roots of P a . Since A'AiA2 = —D, 
we have | Ai A2 1 < 1- So if Ai and A2 are a pair of complex conjugate, then P a i is stable. 
If Ai, A2 are real numbers, by using the inequalities | Ai A2 1 < 1 and P a '(— 1) < < P a '(l) 
and the fact that P a i is a cubic polynomial with positive leading coefficient, we have that 
Ai, A2 lie in the interval (—1, 1) and so P a i is stable. □ 

Proof of Proposition^ From the proof of Proposition^ we have known that —Po(1),Pq(—1) 
and do := —D 2 + AD — B + 1 are the only bifurcation values of the family Pp. Applying 
Proposition [3] to the polynomial, we get that Pp is stable if and only if 

\D\ < 1 and max{-P^(l), P /3 (-l)} < < a(P p ). (10) 

It is easy to see that -Pp{l) = --Po(l) - P, ^(-1) = ^b(-l) - P and a(P p ) = a - (5. 
Hence, (1) implies that Pp is stable if and only if 

\D\ < 1 and max{-P (l), P {-1)} < P < a . (11) 

From (2), we have that the stable interval of Pp exists if and only if 

|D| < 1, d o + Po(l)>0 and a -P (-l)>0. (12) 

By simple computations, we get 

d + P Q (1) = (D + 1)(2 + A — D) and 
do-Po(-l) = (l-D)(2-A + D). 
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Therefore, (3) is equivalent to 

\D\ < 1 and \A-D\<2. 

The equivalence of (a) and (b) is completed. The equivalence of (b) and (c) was actually 
established in the proof of Proposition [3J 

The first statement in item 2 follows immediately from item 1 and conditions (2). 
Finally, we prove the second statement in item 2. By the definition of do, the product 
of two roots of Pa , say Ai and A2, is equal to one. Since Pp is stable for (3 slightly less 
than do, we have that |Ai| = | A2 1 = 1- Because Ai and A2 are different from 1 and —1 (as 
do > max{— P (l), -Po( — 1)}) the result follows. □ 

Acknowledgment. The second and third authors are grateful to Institute of Mathe- 
matics at Academia Sinica in Taipei for hospitality during their visit. 
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Figure 1: The graphs of the bifurcation curves c = ^4^, c = c ^f' and c = ^r, indicated 
simply as Ci , c_x and c, are shown in the (e, c)-plane ("e" is the horizontal axis and 
"c" is the vertical one). The dashed lines in the figures are e = a t+ b i ■ In figures l(i)- 
(m), the stability regions Jt?{F Cfi ) (shaded in black) are shown for three cases: when 
a 2 + b 2 — 1 = 0, > and < 0, all together with b < 0, namely (z) a = 0.6 and b = —0.8, (ii) 
a = 0.2 and b = —1.4, and {Hi) a = 0.1 and b = —0.8. Figure l(iv) corresponds to the 
subcase of (ii) when the stability region has two sides iM" lies inside the strip between 
the dashed lines); here a = 0.85 and b = — 1. 


















A N' 


V 


c, / 







(m) (iv) 
Let us give some remarks on Figure 1. The dashed lines there are e = ± a2 ^ b2 , which 
corresponds to the cases when D — ±1 in Propositions 1 and 2. So the stability region 
Jtr{F Ct e) must belong to the strip between the dashed lines. It is easy to see that the 
intersection of the bifurcation curves c = and c = consists of either one or two 
points depending on whether a 2 + b 2 — 1 is zero or not. In the former case, the point 
of intersection has coordinates (e, c) = (^4^,0), which corresponds to the eigenvalues 
Ai = 1 and A2 3 = aHy/l — a 2 . In the latter case, the two intersection points are M'(e', c') 

and M (e , c ), where c = — 6 , e = ^32, c = — — \ , and 

e" = J ,( ^ b i_ l , which corresponds to the eigenvalues \[ = 1, |A' 2 | = [A3 1 = 1 and A" = \% = 
1, A3 G 1R. Note that the coordinates of M' and M" satisfy the equalities D = 1 and 
D = A + 2 respectively (see conditions in Proposition [31 item 1(a)). Furthermore, the 
condition (5) on existence of stability region Jtv{F Cfi ) implies that e" < e' if a 2 + b 2 — 1 > 0. 
Similar geometric interpretation of the conditions of Proposition [Hand Theorem can be 
done in terms of intersections, say N' and N", of the lines c = c and c = C-\. Note that 
Jtr{F c , e ) has either two "sides" or three "sides" depending on whether both points M" 
and iV" lie outside the strip |e| < a - 2 + b 2 or one of M" and iV" lies inside (one can easily 
see that the points M" and A^" cannot lie simultaneously inside the strip). 
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kc\(e)+c-\(e) 
fe(fe-l) ' 



indicated simply as c, C\ and c_i, are shown in the (e, c)-plane ("e" is the horizontal axis 
and "c" is the vertical one). The dashed lines in the figures are e = a t+ b 2 ■ In figures 
l(i)-(m), the stability regions J non {F c ^ (shaded in black) and Jt T (F^ e ) (shaded in gray) 
are shown for three cases: when a 2 + b 2 — 1 = 0, > and < 0, all together with b < 0, 
namely (i) a = 0.6 and b = —0.8, (ii) a = 0.2 and b = —1.4, and {Hi) a = 0.1 and 
b= -0.8. 
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Figure 3: The graphs of the bifurcation curves c = ~~ ' > c 
fcc -^fc^i)^ 6 ' 1 ; indicated simply as c, c_i and ci, are shown in the (e, c)-plane ("e" is the 
horizontal axis and "c" is the vertical one). The dashed lines in the figures are e = J~^ b 2 ■ 
In figures the stability regions J sym (-Pc,e) (shaded in black), J non (F Cje ) (shaded 

in gray) and Jt T (F Cje ) (shaded in light gray) are shown for three cases: when a 2 + b 2 — 1 = 
0, > and < 0, all together with b < 0, namely (i) a = 0.6 and b = —0.8, (ii) a = 0.2 and 
b = -1.4, and (Hi) a = 0.1 and b = -0.8. 
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Figure 8: Existence of strange attractor. 
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